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Abstract 



The absorption of large bipolarons is investigated using the path-integral 
method. The response of a bipolaron to an external electromagnetic field is 
derived in the framework of the memory-function approach. The bipolaron 
optical absorption spectrum consists of a series of relatively narrow peaks. 
The peculiarities of the bipolaron optical absorption as a function of the 
frequency of the electromagnetic field may be attributed to the transitions 
involving relaxed excited states and scattering states of a bipolaron. 



I. INTRODUCTION 

The optical and kinetic properties of polar and ionic solids are substantially influenced 
by the polaron coupling. Large polarons have been most clearly manifested by investigations 
of the transport phenomena in a magnetic field (see the recent review on polarons0). 

When two electrons (or two holes) interact with each other simultaneously through the 
Coulomb force and via the electron-phonon interaction, either two independent polarons can 
occur or a bound state of two polarons — the bipolaron — can arise (see Refs. on large 
bipolarons and a comprehensive reviewH concerning small bipolarons). Whether bipolarons 
are stable or not, depends on the competition between the repulsive forces (the Coulomb 
interaction) and the attractive forces (mediated through the electron-phonon interaction!. 
Verbist et a/.ii analyzed the large bipolaron using the Feynman path-integral formalismi'0. 
They introduced a "phase diagram" for the polaron — bipolaron system on the basis of a 
generalization of Feynman's trial action and showed that the Pekar-Frohlich coupling con- 
stant as high as 6.8 is needed in 3D to allow for bipolaron formation. Furthermore, in Refs. 

it was shown that the large bipolaron formation is facilitated in 2D as compared to 3D. 
Experimental evidences for bipolarons, e. g. from the data on magnetization and electric 
conductivity in Ti407, as well as in Nao.3V205 and polyacetylene. was discussed by MottS. 

In the framework of the renewed interest in bipolaron theoryBnaflOj a preliminary analysis 
of the absorption of large bipolarons without a magnetic field was givenEl. A variational 
treatment of both spin-singlet and spin-triplet states of large bipolarons (in 2D) and for 
sufficiently strong coupling in high magnetic fields has been presented in Ref. |13|. In Ref. 

the path-integral approach for a bipolaroni was generalized to the case of a bipolaron in 
a magnetic field for all values of the Pekar-Frohlich coupling constant and for all magnetic 



field strengths. It was demonstrated,EJ that the magnetic field favors bipolaron formation. 
The first investigations of the optical absorption spectrum of large bipolarons in a magnetic 
field were performed in Refs. |15|,|16|. 



II. APPROACH 



We investigate here the optical properties of a bipolaron using the path-integral memory- 
function technique, developed in Refs. (see also Ref. |16D. The path- integral variational 
principle gives excellent results for the free energy of the electron-phonon systems (see, e. 
g., Refs. p!9| , pOD . Hence, these models adequately describe the physical properties (including 
the energy spectra) of the polaron and of the bipolaron, respectively. 

The Hamiltonian of two electrons, interacting with longitudinal optical (LO) phonons 
and between each other in an external electric field E {t) , is 



i=i,2 >- 



+ 



£<xj I"! — I"2 



j = l,2 q 



(1) 



with the interaction amplitudesta 



'2^/2' 



1/2 



Txa 



iq 



V 



{- 



h \ 



1/4 



(2) 



where a is the electron-phonon coupling constant, m is the electron band mass, V is the 
volume of the crystal, ^oo is the high-frequency dielectric constant, cjlo is the LO-phonon 
frequency. 

We consider the equation of motion for the vector function R (t) which has the sense of 
the average coordinate response per one electron. 
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Here, the symbol of averaging denotes the path-integral average over trajectories of the 
electrons (cf. Ref. ^ 
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with the action functional 

S [f (t) , f ' it)] = [r it)] - [f (t)] - m [r it) , r' (t)] , (5) 
where Se [ri (t) , r2 (t)] is the action of two interacting electrons in an external electric field. 
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and $ [r (t) , r' (t)] is the "influence phase" of the phonon subsyste 
$[r(t),f'(t)] =-E^/^^/ k (t) - (t)] 
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This influence phase describes both a retarded interaction between different electrons and a 
retarded self-interaction of each electron due to the elimination of the phonon coordinates. 
It contains the Fourier components of the electron density operator 
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and the phonon Green's function 
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T is the temperature. 

We have calculated the optical absorption coefficient for a bipolaron in the memory- 
function approach 
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where c is the velocity of light, n is the refractive index of the crystal, no is the electron 
density. The memory function T [iS) has the formS: 



(e^-* - 1) Im [T^^ (t) (pq it) (0))^J dt. 



(11) 



It is expressed in terms of the two-point correlation function (pq (t) p_q (0))^^ of the electron 
density operators. The correlation function is calculated as the path-integral average [cf. 
Eq. (^] with the model action functional S'o [r (t) ,ff (t) , r' (t) , f'j (t)] : 



/ Df (t) J Df' (t) J Dfj (t) J Df) (t) F exp { (i/h) So [f (t) , (t) , r' (t) , f'^ (t)] } 
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(12) 

The model system consists of two electrons harmonically interacting with two fictitious 
particles and between each other. The model action functional Sq [f (t) ,ff (t) , f' (t) , f'j (t)] 
has the form 
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where M is the mass of a fictitious particle, k, k' and K are the elastic constants. The 
oscillator potentials in the Lagrangian (|1^ imitate the electron-phonon interaction and 
the electron-electron Coulomb repulsion. In Refs. the aforesaid model has been first 
introduced in order to calculate the bipolaron ground-state energy by the Jensen-Feynman 
variational methodi'0. The variational functional of Refs. for the bipolaron free energy 
contains four variational parameters: M, k, k' and K. It is convenient to use instead of 
them four other independent parameters: (i) three eigenfrequencies of the internal bipolaron 
motion fij (z = 1, 2, 3) (see Ref. ||), 
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and (ii) the frequency 
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which is analogous to the Feynman parameter w in the single-polaron problemS. It is seen 
from Eqs. (^) to (p!7[), that the inequality f2i > > v > ^3 is fulfilled for the variational 
frequency parameters. By the variational procedure of Ref. the optimal values of those 
variational parameters are found for the physical system of two electrons interacting with 
the phonon field and between each other. 

In the zero-temperature limit T = 0, we have derived from Eq. (|11]) the analytical 
expressions for the real and imaginary parts of the memory function for a bipolaron in three 
dimensions: 
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From here on, we use the system of units where h = 1, m = 1, cjlo = 1- The factors 
Si {i = 1,2,3) in Eqs. ([ISD , (|l^) have the sense of oscillator strengths corresponding to 
the eigenmodes of the internal motion of a bipolaron. The aforesaid factors are expressed 
through variational frequency parameters of the bipolaron model functional (|15D to ([17D , 

The functions /„ (u) and gn (uj) are given by the expressions: 
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The optical conductivity for a 2D bipolaron is related to that for a 3D bipolaron, which 
is given by the formulas (|18]) and ([19|) , by the scaling relation (cf. Ref. |I]): 



Re cr2D(^, a) 



(24) 



It follows from Eqs. ( p^q ) and (|T9D that the eigenfrequencies 1^2 and appear in the 
optical absorption spectra only in such combinations {k^2 + ^^3) that {k + /) is an even 
integer. This selection rule is determined by the symmetry of these eigenmodes (a schematic 
picture of the internal motion of the bipolaron model system see, e. g., in Ref. |16D. The 
only normal coordinates of eigenmodes with the frequencies ^2 and (^3 (let us denote these 
coordinates as vectors Q2 and Qj) are antisymmetric with respect to the permutation of 
electrons r\ 4=^ r2. Both the exact and model Lagrangians are symmetric with respect to 
this permutation. As a result of this symmetry, we obtain the selection rule 
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(Pq (t) P-q i^))so ^'i- (TO )- The operator of an oscillator coordinate Qjk describes quan- 
tum transitions between energy levels of the corresponding oscillator with the frequency 
Qj (j = 1, 2, 3). Consequently, the aforesaid combinations of coordinates provide the transi- 
tions, which change the energy of the model system by values (^^2 + l^^s) with even [k + I) . 
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III. OPTICAL ABSORPTION SPECTRA 



In order to give a physical interpretation of the calculated bipolaron optical absorption 
spectra, we refer first to a description of those for a polaron in Refs. 0,0. In the polaron 
optical absorption spectra, for intermediate and large a, there is an intense (zero-phonon) 
absorption peak corresponding to a transition from the ground state to the first relaxed 
excited state (RES). The relaxed excited state is created if the electron in the polaron is 
excited while the lattice readapts to a new electronic configurationlii. A shoulder at the 
low-frequency side of the RES peak is attributed to one-phonon transitions through the 
scattering states (ScS): excitations of the polaron system characterized by the presence 
of a finite number of real phonons along with the polaronlll. A broad peak, positioned 
at a higher frequency than the RES peak, is attributed to transitions to Franck-Condon 
(PC) states: internal excited polaron states for which the lattice polarization is that of the 
polaron in its ground stateEl. This picture is well-founded physically and is in agreement 
with the prediction of the aforesaid peaks, which was formulated within the strong-coupling 
approached. 

In Fig. 1, the real ReT (a;) and imaginary ImT {u) parts of the memory function for a 
3D (2D) bipolaron are plotted as a function of frequency for a = 7 (2.971). These values 
of a are close to the minimal possible a for the bipolaron formation: amin = 6.8 (2.9) (see 
Ref. 0). The peaks of the imaginary part of the memory function are positioned near the 
points (explicitly indicated in Fig. 1) 

ujnki = nQi + + + 1, (26) 

where n, k, I are the non- negative integers {k + I = even integer). 

Analogously to the case of polaron optical absorptionEl, the peaks of the bipolaron optical 
absorption spectra (Fig. 2) correspond to: (i) peaks of ImT (a;) and (ii) roots of the equation 
— ReT (lj) = (the crossing points of solid and dotted curves in the panel "a" of Fig. 1). 



Following the physical interpretation developed in Refs. p7| , p!8| , we suggest that the peaks 
corresponding to zeros of [u^ — ReT (cj)] can be attributed to transitions to bipolaron RES. 
Due to a larger number of the internal degrees of freedom for a bipolaron, than that for a 
polaron, there are several types of RES for a bipolaron. 

The positions of peaks of the optical absorption coefficient, corresponding to those of 
ImT (u) , are close to the values Unki determined by Eq. (p6D, so that Unki have the sense of 
the frequencies of transitions from the ground state to certain internal states of a bipolaron. 
Following the physical interpretation developed in Refs. [T7| , [T8| , we can attribute these peaks 
to transitions into Franck-Condon (FC) bipolaron states (in the case when at least one of 
the numbers [n, k, I) ^ 0, i. e., for excited states of a bipolaron). Within the same picture, 
we suggest that the low-intensity peak at = 1 is provided by the transition into the lowest 
scattering state of a bipolaron. It is worth mentioning, that the main resonances in oscillator 
strength are RES. 

In order to classify peaks of the bipolaron optical absorption spectra, we have calculated 
transition frequencies to several bipolaron RES and FC states within the strong-coupling 
(adiabatic) approach. The results are shown in Table 1. The transitions are classified with 
respect to the symmetry of states of the internal bipolaron motion and of the motion of 
a bipolaron as a whole (labelled, respectively, by small and capital letters within standard 
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denotations of states with different orbital moments). "S"' and "s" denote the states with 
the orbital momentum / = 0, while "P" and "p" are those with / = 1. The trial wave 
functions are of the same type as those used for the treatment of polarons in Ref. |18. 



Table 1. Transition frequencies (in units of the LO phonon frequency) to 
bipolaron RES and FC states calculated within the adiabatic strong-coupling 

theory. 



a 


RES 


FC 




pS 


sP 


pP 


pS 


sP 


pP 
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2.2171 


5.9877 


5.7282 


5.4633 


9.9396 


12.5210 


8 


2.8957 


7.8207 


7.4818 


7.1328 


13.0008 


16.3627 
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3.6648 


9.8981 


9.4686 


8.9364 


16.2629 


20.4673 


10 


4.5243 


12.220 


11.6895 


11.1019 


20.1550 


25.4065 



Comparing the peak positions of Fig. 2 with the results of Table 1 (which are obtained in 
the framework of the adiabatic approach, and therefore can be applied only approximately), 
we have attributed several optical absorption peaks to transitions to RES and FC states of 
a definite symmetry, as shown at the figure. It should be mentioned, that every RES peak 
has a lower frequency, than the corresponding FC peak, since the lattice relaxation leads to 
a lowering of the energy of the electron-phonon system. 

The linewidth of peaks of bipolaron optical absorption spectra is determined by the 
bipolaron recoil in the scattering process. Every FC peak is described by a function 
Qn+k+i — ^nki) of Eq. (p!9D , where gn (w) is given by Eq. (|^). Consequently, the linewidth 
of a FC peak can be estimated, using the following characteristic of the function (p2|): 



Jo 9n{uj)duj 

Performing the integrations over uj analytically, we find that 

B f 3^ 



^ \9n (^)] = ^ (^^ + 2 

where the parameters A and B are given by Eq. (H). The ratio | is of the same order as 
[2 (M + 1)] , which is used for estimation of the bipolaron effective massed. Therefore, the 
characteristic linewidth of a FC peak corresponding to a definite frequency uj = uOnki can be 
qualitatively estimated as 

Since the bipolaron effective mass is larger than that of a polaron, the bipolaron optical 
absorption spectrum consists of a series of peaks which are narrower in comparison with 
those of the polaron absorption spectrum. 

Let us define N (cj, Acu) as the number of possible combinations of (n, k, I) for which 
the frequencies {uJnki} occupy the interval uo ^ uOnki ^ ^ + Acu. The higher is uo (at a fixed 
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interval length Au), the larger is N {u, Au) . Hence, the density of the FC peaks increases 
with increasing frequency. Starting from a definite frequency range (where {u, Au) oc 
Tpc i^) Au), FC peaks merge into a continuous band. 

In the limiting case u ^ 1, the imaginary part of the memory function (|19D asymptoti- 
cally behaves as 

ImTHU^oc-^^^/l (28) 

Hence, the bipolaron absorption coefficient in the limiting case of high frequencies has the 
same asymptotic behavior as the polaron absorption coefficients^ (this asymptotic behavior 



was also derived in Ref . for the weak-coupling polaron optical absorption) : 



r(c.)U,cx.;-^/^ (29) 

This behavior has a clear physical explanation: in the high-frequency limit, only one-phonon 
scattering processes give a contribution into the optical absorption by the electron-phonon 
system. In this limit, the optical absorption coefficients for bipolarons and polarons at 
arbitrary coupling strength are described by one and the same asymptotic dependence. It 



was emphasized in Refs. 17,15,24, that the frequency dependence (M) of the polaron optical 



absorption coefficient at u ^ 1 differs from the Drude-like dependence r*^^'^'^'^^^ (to) | oc 

For the bipolaron, the optical absorption spectrum turns into a continuous band at 
substantially higher frequencies, than for the polaron, due to a comparatively large effective 
mass of a bipolaron. This frequency region is not presented in Fig. 2. The bipolaron 
absorption coefficient in this region is very small when compared with that shown in Fig. 2. 



IV. SUM RULES AND CONCLUSIONS 

From the general analytical properties of the memory function, for the real part of the 
optical conductivity per one electron 

Rea (uj) = ^ ^ (30) 

the following sum rule was derived in Ref. |2^: 



where m* is the polaron effective mass. For the bipolaron optical absorption, the same 
equation (|3T|) is valid with m* the bipolaron effective mass. Our numerical test of the sum 
rule ( ^TD has confirmed, that Eq. ( ^T]) is fulfilled within the chosen relative accuracy (up to 

10-3). 



In Ref. 26, the ground state theorem for a polaron has been derived which relates the 



polaron ground state energy to the first moment of the optical absorption spectra. The 



extension of this theorem to the bipolaron can be performed in the same way as in Ref. |26 . 
The analog of Eq. (25) from Ref. ^ for a bipolaron is 
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where ai and 02 are two arbitrary values of a, r/ = Eoo/eo is the ratio of the high-freauency 
and static dielectric constants, Xjj i^^cniV) is the current-current correlation functioned. For 
a bipolaron, we choose 1] and both ai and 02 within the region of the bipolaron stability. 
The function lw.Xjj {^) for a bipolaron is expressed through the memory function T {uj) as 

Imx,-,- = ^. (33) 



In Fig. 3, the ground-state theorem (32) is illustrated numerically. The solid curve shows 
the difference between the ground state energies Eq {a2,ri) — Eq (ai,?7) (calculated by the 
variational method using the variational functional of Ref. ^) for two values of a at 77 = 0. 
The value ai is taken ai = 7. The dashed curve shows the values of the right-hand side of 
Eq. (|32|). As seen from Fig. 3, there exists an excellent agreement between the quantities 
entering the independently calculated left-hand and the right-hand sides of the equation 



In the present paper, we have treated the bipolaron optical absorption within the 
memory-function path-integral approach. The derived optical absorption spectra demon- 
strate a rich structure of relatively narrow peaks, which is clearly related to the internal 
states of a bipolaron. We have attributed those peaks of the bipolaron optical absorption 
spectra to transitions from the ground state to (i) scattering states, (ii) relaxed excited 
states, (iii) Franck-Condon states. Every RES peak is shifted to a lower frequency with 
respect to the corresponding FC peak. This interpretation is performed within a unified 



approach to the polaron and bipolaron problems, which takes its origin in Refs. p!7| , p!8| . The 
sum rules, developed for the polaron optical conductivitySS, are extended to the case of a 
bipolaron. 
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Figure captions 



Fig. 1. Real [panel "a"] and imaginary [panel "b"] parts of the memory function T (u) 
for the optical absorption coefficient of a 3D (2D) bipolaron for a = 7 (2.971), t] = 0.0037. 
The optimal values of the variational frequencies are indicated in the figure. 

Fig. 2. Optical absorption spectra of a 3D (2D) bipolaron for a = 7 (2.971), t] = 0.0037 
[panel "a"], for a = 8 (3.395), r] = 0.023 [panel "b"], and for a = 9 (3.820), r] = 0.023 [panel 
"c"]. 

Fig. 3. Illustration of the ground-state theorem (^) for the bipolaron optical absorption. 
The index "1" denotes the ground state energies calculated by the variational method [the 
left-hand side of Eq. (0)], while the index "2" denotes the ground state energies calculated 
by integration of the bipolaron optical absorption spectra [the right-hand side of Eq. (^)]. 
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Fig. 1. J. T. Devreese, S. N. Klimin, and V. M. Fomin 
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Fig. 2 J. T. Devreese, S. N. Klimin, and V. M. Fomin 




Fig. 3 J. T. Devreese, S. N. Klimin, and V. M. Fomin 



